arXiv:1504.03904v2 [quant-ph] 6 May 2015 


APS 


Experimental characterization of the Gaussian state of squeezed light obtained via 

single-passage through an atomic vapor 

P. Valente, A. Auyuanet, S. Barreiro, H. Failache, and A. Lezamaj^ 

Instituto de Fisica, Facultad de Ingenieria, Universidad de la Republica, 

J. Flerrera y Reissig 565, 11300 Montevideo, Uruguay 
(Dated: May 7, 2015) 

We show that the description of light in terms of Stokes operators in combination with the as¬ 
sumption of Gaussian statistics results in a dramatic simplification of the experimental study of 
fluctuations in the light transmitted through an atomic vapor: no local oscillator is required, the 
detected quadrature is easily selected by a wave-plate angle and the complete noise ellipsis recon¬ 
struction is obtained via matrix diagonalization. We provide empirical support for the assumption 
of Gaussian statistics in quasi-resonant light transmitted through an ®^Rb vapor cell and we illus¬ 
trate the suggested approach by studying the evolution of the fluctuation ellipsis as a function of 
laser detuning. Applying the method to two light beams obtained by parting squeezed light in a 
beamsplitter, we have measured entanglement and quantum Gaussian discord. 
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I. INTRODUCTION 

It is well established that the propagation of linearly 
polarized continuous-wave laser light through a near res¬ 
onant atomic vapor can result in squeezing of the field 
with polarization component orthogonal to the incident 
field polarization (vacuum squeezing) [TH5]. Such phe¬ 
nomenon has been understood as the consequence of the 
nonlinear response of the atomic medium giving rise to 
the effect of polarization self-rotation (PSR) of elliptically 
polarized light. PSR is the consequence of the third order 
susceptibility of the atomic sample. Using this suscepti¬ 
bility, it can be shown that the input-output relations for 
the quantum held operators [2] result in vacuum squeez¬ 
ing of the field polarized in the direction orthogonal to 
the incident field polarization (hereafter referred as “vac¬ 
uum polarization”). 

Initial work on this subject has centred the attention 
only on the squeezed polarization component of the held. 
From the experimental point of view, this required an 
auxiliary local oscillator to perform the balanced homo¬ 
dyne detection of the squeezed light. Usually, the strong 
linear polarization component of the held has been used 
as a local oscillator mis]. 

In a previous article [6], we have adopted a different 
point of view in which the two polarization components 
of the light beam are simultaneously considered. From 
this perspective, the squeezing effect is seen as a modih- 
cation in the fluctuations of the polarization of the light 
i.e. the Stokes operators associated to the held. A re¬ 
duction of the variance of a Stokes operator below the 
corresponding standard quantum limit (SQL) indicates 
polarization squeezing [7]. Experimentally, the main ad¬ 
vantage of using polarization modes for the investigation 
of the quantum properties of light is the elimination of 
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the requirement of a local oscillator. 

Polarization squeezing and entanglement were hrst ob¬ 
tained using two independent optical parametric ampli- 
hers [SE] and a cold atomic cesium cloud inside an op¬ 
tical cavity [ini? Taking advantage of the fact that 

polarization squeezing does not require a shared local 
oscillator, it was recently used for entanglement distri¬ 
bution [12j and long distance quantum communication 

m- 

In this article we go a step forward in the examination 
of the huctuation properties of the polarization state of 
laser light transmitted through an atomic sample by in¬ 
corporating the additional assumption of the Gaussian 
statistics of the quantum state of the light held. These 
states have a very elegant mathematical description and 
are simple to manipulate. Also, Gaussian states of light 
are commonly encountered: the vacuum and the coher¬ 
ent states are Gaussian and so are the squeezed vacuum 
and the displaced squeezed vacua. In addition, a wide 
range of held evolutions are known to preserve the Gaus¬ 
sian nature of the held state [14] . In particular, the PSR 
effect results in squeezed vacuum of the vacuum polariza¬ 
tion component. Based on the above, we have made the 
working assumption of the Gaussian nature of the light 
state. We provide experimental evidence in support of 
this assumption. 

The purpose of this article is to show that by combining 
the description of the held in terms of Stokes operators 
and the hypothesis of Gaussian statistics one can easily 
have access to full information about the light state. We 
illustrate this possibility by studying laser frequency de¬ 
pendence of the squeezed and anti-squeezed quadrature 
variances and the corresponding quadrature angles. In 
the last part of the paper, we apply Gaussian analysis of 
Stokes operators measurements to the determination of 
entanglement quantihers and quantum Gaussian discord 
on two light beams separated in a beamsplitter. 
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II. BACKGROUND 
A. Stokes operators 

A quasi-monochromatic light beam propagating in the 
z direction is associated to the annihilation operators 
Ux and Uy corresponding to the x and y linear po¬ 
larization components respectively. Alternatively, one 
can refer to a different polarization basis and introduce 
ai( 2 ) = [ax + [—)ay)/'/2 corresponding to the polariza¬ 
tion components along the main diagonals of the cc, y 
plane or a± = (ax ± iay)f'/2 corresponding to the two 
circular polarizations components. 

In order to describe the polarization state of the held 
it is convenient to introduce the Stokes operators: 


al,ax + a/yay = nx+ny 

(la) 

aiax - alay = Ux - riy 

(lb) 

al,ay + a^ax = ni — n 2 

(Ic) 

i [a^yax — a\ay) = n_ — 

(Id) 


where rii, = ala^, is the photon number operator associ¬ 
ated to a given polarization. The Stokes operators S 1-3 
obey the angular-momentum-like commutation rule: 

[Si,Sj] = 2ieijkSk, ( 2 ) 

where eijk is ±1 depending on the parity of the i,j,k 
permutation. 

It is convenient to introduce a generalized Stokes op¬ 
erator in the 2 — 3 plane parametrized by the angle 0: 

Se = ne+ — ng^ = S2cos[0) + S3sin[9) (3) 

where ng± are the photon number operators correspond¬ 
ing to the held operators: ag± = [ax ± e^^ay)fy /2 which 
correspond to two orthogonal elliptical polarization com¬ 
ponents. It can be easily shown that Sg^o = S 2 and 
Sg=Tr /2 = S 3 . The generalized Stokes operators verify 
the commutation rule: 

[Sg,Sg>] = 2iSisin[0' -0) (4) 

In this article we are concerned with a held that is in¬ 
tense and essentially polarized along the x direction. Let 
d = [ax) be the expectation value of the held operator 
for the X polarization. We assume that |q!| 3> | (a^) | ~ 0 
and that the huctuations of the x polarization compo¬ 
nents of the held are small compared to a. Under such 
conditions the Stokes operators can be approximated as: 


^0 ^ lap (5a) 

S '2 ~ a*ay -I- aa'^y (5b) 

S 3 ~ i [aa^y — a*ay) (5c) 

Sg ~ a*aye~^^ + aa^yC''^ (5d) 


with 

\Sg^ Sgi] 2i\a\'^sin[0'— 9 ) ( 6 ) 

In particular, 

[Sg,Sg+^^]~2l\a\^ (7) 

which leads to the uncertainty relation: 

A5,A5(e+.) > lap (8) 

Introducing Xg = Sg/\a\ and Pg = Sg+^/\a\ we have: 

[Xg,Pg] = 2l (9) 

which is the usual commutation rule for orthogonal 
quadratures. In consequence these operators must obey 
the Heisenberg uncertainty relation: 

AXgAPg ^ 1 (10) 

The observables Sg and Sg+^ as well as |ap and there¬ 

fore Xg and Pg, are accessible to measurement. In this 
paper we consider a phase-space description of the held 
state in the X, P (or Sg, Sg+^) plane. 



FIG. 1: (Color on-line) Experimental scheme for the measure¬ 
ment of the Stokes operator Sg. The angle 9 is selected by 
the orientation of the half-wave plate axis. PBS: polarization 
beam-splitter. PD: photodetector. 

The measurement of the generalized operator Sg can 
be achieved using the arrangement shown in Fig. A 
hxed quarter-wave plate positioned with its principal axis 
oriented at 7r/4 with respect to the axis x is followed 
by a half-wave plate whose fast axis forms an angle of 
tt/S -I- 9/4: with respect to x. With this arrangement, 
the two photodectectors outputs are respectively propor¬ 
tional to n6/_|_ and ng_ i.e. the photon numbers corre¬ 
sponding to two orthogonal elliptical polarizations. In 
the case of our experiment, the axis of the two ellipsis 
are oriented at ±7r/4 with respect to the incident po¬ 
larization. In consequence, the mean value of the two 
photodetectors outputs are balanced. This arrangement 
is simpler than the one that we have previously used [6] 
which required a tilted waveplate and the calibration of 
the phase difference between orthogonal linear polariza¬ 
tion as a function of the tilt angle. 
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B. Gaussian states 

In recent years, quantum information processing with 
continuous variables has received large attention M- A 
fundamental role is played by quantum systems prepared 
in Gaussian states [T^IITC] . Such states have been used to 
demonstrate quantum correlations between two separate 
systems. Over the past decade, a significant amount of 
work has been devoted to the demonstration and mea¬ 
surement of entanglement using continuous variables [ni- 
[27]. More recently [JSHSl], Gaussian states were used in 
the determination of the quantum discord between two 
systems. 

A quantum state describing a continuous variable is 
said to be Gaussian if its Wigner function defined over 
the X, P plane is Gaussian. A two-dimensional Gaussian 
function is entirely defined by the covariance matrix: 

(AX^ AXP\ 

[axP AP2 J 

where for any two operators Y, Z we use AZYZ = 

i(FZ + zr)-(y)(z). 

For a given 0, it is easily shown that the covariance 
matrix is given by: 


/ AXl 


\\AXl^^^-AXl_ 
AX- 



where \ota\’^ and \ab\- are given by the shot noise levels 
on modes a and b respectively. 

Using the setup shown in Fig. the variances 

A{Sga ± Sg'b)-', Corresponding to the sum and 
difference of the two balanced detectors outputs, 
can be measured. With the help of the identity, 
ASeaSg.b = [AiSea + Sg^b? - AiSea - Se'b f]/4: and 
the previously determined values of \aa\ and lahl the 
two-mode coefficients are readily computed. 

Williamson’s theorem [32] ensures that the covariance 
matrix can be diagonalized through the application of a 
symplectic transformation. Also, the covariance matrix is 
characterized by four invariants under symplectic trans¬ 
formations. The two eigenvalues of the diagonal form 
are the (degenerate) symplectic eigenvalues and v- 
which can be expressed in terms of the symplectic in¬ 
variants [33| . The following properties derive from the 
symplectic invariants and eigenvalues. 


1. Consistency 


A covariance matrix must satisfy some constrains to 
represent a physical Gaussian state. In order to ver¬ 
ify the Heisenberg uncertainty relation, the symplectic 
eigenvalues Vk must verify [53] : 


> 1 


(15) 


All terms in (121 can be experimentally measured. 


Once the covariance matrix is known for a given choice 
of the angle 0 it can be easily calculated for different 
angles through a rotation transform (see Eq. [^. As 
the value of 9 is varied, the variance AXg describes 
an ellipsis. For the proper choice of 0 corresponding 
to the orientation of one of the ellipsis main axis, the 
covariance matrix becomes diagonal. If AXg < 1 (see 
Eq. 101 the field is squeezed and the angle 9 identifies 


the corresponding quadrature. 


For two modes (a and h) a Gaussian state is entirely 
characterized by the covariance matrix: 


AXl AXaPa 
APaXa AP- 
AXbXa AXbPa 

yAPbXa APbPa 


AXaXb AXaPb\ 

APaXb APaPb 
AX- AXbPb 
APbXb AP- / 


(13) 


All the coefficients in (131 can be experimentally mea¬ 
sured using the setup shown in Fig. The coefficients 
referring to a single mode are determined as in (121. The 


two-mode coefficients of the covariance matrix are of the 
form: 

ASgaSgib 


AYaZb = 




(14) 


2. Entanglement 

The positivity of the partial transposed matrix is a nec¬ 
essary condition for separability [3S] . It is also a sufficient 
condition in the case of two-mode Gaussian states [36] . 
If and £'+ with (£'-|_ > f'-) are the symplectic eigen¬ 
values of the partially transposed covariance matrix, the 
state is entangled if :>_ < 1. In addition, the amount of 
entanglement can be quantified through the logarithmic 
negativity LN = — ^ogivk) (the sum extending over 
all Vk < 1) or LN = 0 if f'fe > 1,VA: [57] . 


3. Quantum Gaussian discord 

Entanglement does not account for all quantum corre¬ 
lations. In fact, quantum correlations can exist for sep¬ 
arable mixed states. The quantum discord |38j intends 
to quantify quantum correlations beyond entanglement. 
Initially defined for discrete quantum systems, its def¬ 
inition was extended to Gaussian states [55] SO]. The 
optimality of the quantum Gaussian discord (QGD) is 
discussed in [4T]. The QGD can be directly computed 
from the four symplectic invariants of the covariance ma¬ 
trix [30] . 
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III. EXPERIMENT 

A scheme of the experimental setup is presented in Fig. 
m The details of the laser source and the atomic cell 
environment have been previously described in [^. We 
remind the essential features. A 795 nm laser beam with 
approximately 40 mW of total power is linearly polarized 
and focussed with a 50 mm focal length at the center of 
a 5 cm long uncoated glass vapor cell containing isotopi- 
cally pure ®^Rb surrounded by a magnetic shield. The 
cell is heated to 50° C. The laser beam is re-collimated 
after the cell. The laser frequency can be scanned around 
the Rb D1 transitions. 



FIG. 2: (Color on-line) Experimental setup used for the char¬ 
acterization of a two-mode state corresponding to the out¬ 
put of the ~ 50% non-polarizing beam-splitter (BS). On each 
output a Stokes operator measurement setup is implemented. 
The balanced detector outputs can be either added or sub¬ 
tracted for global operator mesurements. 

In the first set of experiments described below the 
beam splitter (BS) shown in Fig. was not present 
thus directly realizing the detection arrangement pre¬ 
sented in Fig. [2 To determine the shot-noise level, a 
linear polarizer is introduced on the beam path before 
the wave-plates. As a consequence, the light polariza¬ 
tion perpendicular to that of the incident field is blocked 
and replaced by vacuum. In this case, the balanced de¬ 
tector output variance corresponds to ASg = lap which 
represents the shot-noise level for the fluctuations of the 
Stokes operators Sg. In subsequent measurements, with 
the polarizer removed, the shot-noise level is scaled in 
proportion to the light intensity to account for the small 
attenuation introduced by the polarizer. 

In a second set of experiments the non-polarizing 
beamsplitter BS was introduced. An approximately 50% 
transmission and reflection BS was used. Two separate 
balance detections were implemented to measure Stokes 
operators on both output ports of the beamsplitter. As 
a consequence of the BS being non-polarizing, the in¬ 
tense polarization component of the field is the same on 
the two outputs. The wave-plate angles are referred to 
this common direction. The AC outputs of the two bal¬ 


anced detectors are either added or subtracted to access 
the fluctuations of global observables of the system. The 
total signal is sent to an oscilloscope for a time domain 
study of the fluctuations or to a spectrum analyser (SA) 
that records the noise power within a given bandwidth 
around a selected central frequency. The shot-noise levels 
corresponding the each of the light beams are indepen¬ 
dently recorded. 

For a given combination of the waveplates, the laser is 
scanned around the D1 transitions of ®^Rb and the cor¬ 
responding fluctuations data is recorded. Using the data 
for the four choices of 6* = —7r/4, 0,7r/4,7r/2 the covari¬ 
ance matrix (Eq. Ill is determined as a function of the 
laser detuning. For each detuning the covariance matrix 
can be brought into its diagonal form and the minimum 
and maximum noise variances computed as well as the 
corresponding quadrature angles. 


IV. RESULTS 
A. Gaussianity test 

We have initially tested our assumption of Gaussian 
statistics of the light fluctuations. To assert the Gaus¬ 
sian nature of the field state is not a simple task [15] . We 
have approached this question empirically by examining 
the data corresponding to the balanced detector output 
as a function of time for different choices of the Stokes op¬ 
erator and the field. For each Sg, we have recorded 40 os¬ 
cilloscope traces consisting of 2500 data points taken over 
a 250 pLS time interval. Since the temporal record of the 
fluctuations is plagued by technical noise excess, specially 
at low frequency, we have numerically band-pass filtered 
the data string to extract the fluctuations corresponding 
to a specific frequency interval [fl — A/2, U -|- A/2] where 
U is the noise central frequency and A the bandwidth. 
The filtering procedure consisted in Fourier transform¬ 
ing the data string, replacing by zeros the points corre¬ 
sponding to frequencies whose absolute value lies outside 
the frequency interval and performing the inverse Fourier 
transform. 

The histogram of the filtered time domain fluctuations 
are presented in Fig. [^for different choices of the Stokes 
operator Sg. The histogram of the corresponding shot 
noise fluctuations is also shown. All histograms are well 
fitted by a normalized Gaussian function with the cor¬ 
responding variance. The data in Fig. [^corresponds to 
a maximum squeezing of 1.6 dB with no correction for 
losses. 

An additional test of the Gaussian statistics of the 
Stokes operators fluctuations results from the computa¬ 
tion of the stochastic moments. For a Gaussian distribu¬ 
tion, the odd moments are all zero and for p even they 
obey: (x^) = (p — 1)!!. Figured shows the com¬ 

parison of the statistical moments of the recorded data 
to the corresponding value for a Gaussian distribution 
for moments up to p = 6. The odd moments are zero 
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FIG. 3: (Color on-line) Symbols: Normalized histograms of 
the measurement of the filtered (fl = 3 MHz, A = 400 
kHz) fluctuations of the Stokes operators. Solid lines: Gaus¬ 
sian functions with the same variance than the corresponding 
data. Dashed line: Residual electronic-noise level. The laser 
is tuned near the ®^Rb Fg = 2 ^ Fe = 2 transition at the 
position of maximum squeezing. 


within the measurement uncertainty and the even mo¬ 
ments differ from the expecting Gaussian value in less 
than 8 %. 
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FIG. 4: (Color on-line) Normalized statistical moments for 
the measurements of Stokes operators. The normalized mo¬ 
ment of order p is defined as: (a:^)/ (p — l)!!j — Sp 

where £p = 1 if p is even and zero otherwise. Here x repre¬ 
sents the result of a measurement. The normalized moments 
are zero for a Gaussian distribution. 

Supported on these observations, we conclude that 
the assumption of Gaussian nature of the field state is 
reasonable in the context of this work. 
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B. Single Stokes operator 

Figure presents the noise power (recorded with a 
spectrum analyser) around 2.7 MHz (resolution band¬ 
width 100 KHz) of a single Stokes operator as a func¬ 
tion of laser detuning. In Fig. Se was chosen to cor¬ 
respond to the maximum squeezing occurring near the 
F = 2 —>■ F' = 2 transition. Smaller squeezing is also 
observed around the F = 1 ^ F' transitions. For other 


laser detunings the fluctuations of Sg are above the shot- 
noise level. It is worth noticing the broad frequency range 
for which the light fluctuations are affected by the atomic 
transitions. It extends over a range of around 17 GHz, 
much larger than the total D1 absorption spectrum of 
rubidium. A slow decrease of the excess noise structure 
with the laser detuning from atomic resonances was pre¬ 
dicted by the numerical simulations reported in | 3 ]. 



-5 0 5 10 

Laser detuning (GHz) 


FIG. 5: (Golor on-line) a) Reference-cell saturated-absorption 
signal used for frequency calibration. The hyperfine transi¬ 
tions of the ®^Rb D1 line are indicated, b) Gircles: Noise 
power at 2.7 MHz for the Stokes operator Sg presenting max¬ 
imum squeezing around the F = 2 ^ F' = 2 transition in a 
cell containing isotopically pure ®^Rb. Triangles: Shot noise 
level. 


C. Phase-space description of one field 


In order to have access to the covariance matrix, 
we have recorded with a spectrum analyser (cen¬ 
tral frequency 2.7 MHz, resolution bandwidth 100 
kHz) the fluctuations of the four Stokes operators 
Sg- 7 r/ 4 ,-Sg, S'g_|_.,r /45 >S'e-i- 7 r/ 2 - For Convenience we chose Sg 
to correspond to S' 3 . The covariance matrix was con¬ 
structed according to ( 12 ) and diagonalized to identify 


the minimum and maximum noise variances and the cor¬ 
responding quadrature angle. The results are presented 
in Fig. |^(see also |43] ). Similar spectra of the minimum 
and maximum noise variances were presented in [Si]. 
Notice the considerable simplification arising from the 
combination of Stokes operators detection and Gaussian 
analysis. In 011], a sophisticated active stabilization of 
the phase of the local oscillator was implemented in order 
to track the field quadrature corresponding to maximum 
or minimum noise. In |S] this was achieved manually for 
a discrete number of laser detunings. Also these articles 
did not report the determination of the quadrature an¬ 
gle corresponding to squeezing possibly due to difficulties 
in the local oscillator phase calibration in the interfero¬ 
metric setup. By contrast, in our setup there are no 
interferometric instabilities affecting the choice of Sg , no 
stabilization loop is required since the data is recorded 
for a fixed value of the angle 9 and the determination 
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of the noise variances and corresponding quadrature an¬ 
gles readily arise from the numerical covariance matrix 
diagonalization. 



I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I 

-4 -2 0 2 4 


Laser detuning (GHz) 

FIG. 6: (Color on-line) (a) Quadrature noise power normal¬ 
ized to the shot-noise level as a function of laser detuning. 
Squares: Minimum quadrature noise. Circles: Maximum 
quadrature noise. Solid: Reference saturated absorption sig¬ 
nal with relevant ®^Rb D1 transitions indicated, (b) Quadra¬ 
ture angle corresponding to the minimum noise. (See also 

m)- 


The results presented in Fig. [^reveal interesting fea¬ 
tures about the noise structure. As previously observed 
EHS], squeezing is present around the F = 2 ^ F' = 1 
and F = 2 —>■ F' = 2 albeit for different quadrature an¬ 
gles. Interestingly enough, for the laser detunings where 
squeezing does not occur, the minimum noise variance 
remains close to the shot-noise level. On the other hand 
the maximum noise variance presents excess noise above 
the limit imposed by the Heisenberg uncertainty relation. 

Perhaps the more interesting feature in Fig. is the 
existence of a laser frequency between the two hyperfine 
transitions for which the maximum and minimum noise 
variances are nearly equal and very close to the shot- 
noise level. The existence of this singular point in the 
noise spectrum can be related to the cancellation of the 
PSR effect around this position [1]. It is interesting to 
notice that not only the squeezing is suppressed but very 
little excess noise is added by the atomic interaction. The 
squeezing occurs on either side of such noise cancellation 
point. 

Our method readily allows the observation of the ro¬ 
tation of the Stokes operator quadrature noise ellipsis as 
a function of laser detuning |43) . A steep variation of 
the orientation of the noise ellipsis over 72 degrees oc¬ 
curs around the noise cancellation point. It was recently 
suggested [HI SS] that the control of the orientation of 
squeezed light noise ellipsis could be applied to precision 
improvement in gravitational wave detection interferom¬ 
etry. Rotation of the noise ellipsis was recently observed 
in squeezed light generated via four-wave-mixing in an 
atomic system |46j . 


D. Two-field state characterization 



Laser detuning (GHz) 


FIG. 7: (Color on-line) Two-mode state quantifiers for the 
two beamsplitter outputs as a function of laser detuning. Cir¬ 
cles: minimum partial transpose eigenvalue {MPTE = V-). 
Squares: logarithmic negativity (LN). Triangles: quantum 
Gaussian discord (QGD). 


We have recorded the covariance matrix of the com¬ 
pound system for each value of the laser detuning follow¬ 
ing the procedure described in Section jlT^ 

Since the experimentally determined covariance matrix 
is affected by noise, we have tested each covariance ma¬ 
trix for consistency and physical significance. We have 
discarded all the covariance matrices that did not satisfy 
the physical consistency requirement (Eq. 15). The data 
complying to this requirement was used for symplectic 
invariants determination and computation of the sym¬ 
plectic eigenvalue = f'_, the logarithm negativity and 
the quantum Gaussian discord between modes a and b 
conditioned to measurements on b [40] . 

The results are presented in Fig. [^ In spite of the 
significant amount of noise, the positivity of the par¬ 
tial transpose criterion and the logarithmic negativity 
indicate entanglement near the laser frequency range for 
which squeezing is present on single held Stokes opera¬ 
tors (see Fig. [^. Notice that nonzero quantum Gaussian 
discord occur in all the laser frequency range shown in 
Fig. [^indicating that the two-mode state presents some 
amount of nonclassical correlations even when it is sepa¬ 
rable. 

As a hnal remark, it is worth mentioning that the ex¬ 
perimental procedure described above is independent of 
the intensity balance of the two outputs of the BS. Such 
balance will nevertheless have an inhuence on the entan¬ 
glement and quantum Gaussian discord observed. 


V. CONCLUSIONS 


We have shown that the description of the held in 
terms of polarization (Stokes) states in combination with 
the assumption of Gaussian statistics results in consider¬ 
able simplihcation of the study of the quantum state of 
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light transmitted through an atomic sample: The experi¬ 
mental setup is dramatically reduced. All interferometric 
instabilities are removed. The selection of the detected 
phase space quadrature corresponds to the rotation angle 
of a single waveplate. Finally, the complete state (noise 
ellipsis) reconstruction is computed through matrix di- 
agonalization. 

Taking advantage of this simplicity, we have presented, 
for the first time to our knowledge, the complete evo¬ 
lution of the noise ellipsis describing light fluctuations 
around a given central frequency as a function of laser 
detuning. In the case of two light modes, obtained by 
parting a squeezed light beam in a beamsplitter, the use 


of Gaussian states formalism in combination with Stokes 
operator observations, allowed us to demonstrate entan¬ 
glement and measure the quantum Gaussian discord in¬ 
dicating that the generated light states could be suitable 
for quantum information processing. 
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